Abstract. Studies of the Fundamental Plane of early-type galaxies, from small to in- 
termediate redshifts, are generally carried out under the guiding principle that the Fun- 
damental Plane reflects the existence of an underlying mass-luminosity relation for such 
galaxies, in a scenario where galaxies are homologous systems in dynamical equilibrium. 
In this paper we reexamine the question of whether a systematic non-homology could 
be partly responsible for the correlations that define the Fundamental Plane. We start 
by studying a small set of objects characterized by photometric profiles that have been 
, pointed out to deviate significantly from the standard i? 1 / 4 law. For these objects we 
£N) ■ confirm that a generic i? 1 /™ law, with n a free parameter, can provide superior fits (the 
i-C ■ best-fit value of n can be lower than 2.5 or higher than 10), better than those that can be 
Ph ' obtained by a pure i? 1 / 4 law, by an .R^+exponential model, and by other dynamically 
,— ' \ justified self-consistent models. Therefore, strictly speaking, elliptical galaxies should not 
be considered homologous dynamical systems. Still, a case for weak homology, useful for 

i— I ■ 

■ the interpretation of the Fundamental Plane, could be made if the best-fit parameter n, 
00 ' 

' as often reported, correlates with galaxy luminosity L, provided the underlying dynami- 



(N 
(N 



cal structure also follows a systematic trend with luminosity. We demonstrate that this 



, statement may be true even in the presence of significant scatter in the correlation n(L). 

■ Preliminary indications provided by a set of "data points" associated with a sample of 14 

' galaxies suggest that neither the strict homology nor the constant stellar mass-to-light 
1^ solution are a satisfactory explanation of the observed Fundamental Plane. These con- 
+3 . elusions await further extensions and clarifications, because the class of low-luminosity 
CS • early-type galaxies, which contribute significantly to the Fundamental Plane, falls outside 
m ^ ' the simple dynamical framework considered here and because dynamical considerations 

' should be supplemented with other important constraints deriving from the evolution of 

b : 

5t i stellar populations. 
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1. Introduction 

In the three-dimensional space of observed parameters loga , \ogR c , and SB C = 
— 2.51og(/) c (where (I) c = L/27ri? 2 ), early-type galaxies are approximately located on 
a plane, called the Fundamental Plane (hereafter FP; Dressier et al. 1987; Djorgovski & 
Davis 1987; Bender, Burstein, & Faber 1992), represented by the best-fit relation 
logi? = a\og<7 + [3SB c + j . (1) 
The central velocity dispersion <j is often defined as the velocity dispersion referred 
to an aperture radius of R c /8 (J0rgensen, Franx, & Kjaergaard 1993); the quantities 
R c and L are the empirically determined half-luminosity (effective) radius and total 
absolute luminosity. The coefficients a, (3, and 7 depend slightly on the photometric band 
considered. By measuring R e in kpc, cr in km s _1 , and SB C = 42.0521 — 2.5 log(L/27ri? 2 ) 
in mag/ arc sec 2 , where L is expressed in units of the solar Blue luminosity, reported values 
are a = 1.25 ± 0.1, (3 = 0.32 ± 0.03, 7 = -8.895 (J0rgcnsen, Franx, & Kjaergaard 1993; 
Bender ct al. 1998; note that these coefficients, in particular this choice of 7, refer to a 
long cosmological distance scale). The FP is characterized by a small scatter in R c , on 
the order of 15%, remarkably independent of the position on the FP itself. 

Studies of the FP of early-type galaxies, from low to intermediate redshifts (for this 
latter case, see van Dokkum & Franx 1996; Kelson et al. 1997; van Dokkum et al. 1998a,b; 
Bender ct al. 1998; Pahre, de Carvalho, & Djorgovski 1998; Pahre, Djorgovski, & de Car- 
valho 1998; J0rgensen et al. 1999; Treu et al. 1999; Kelson et al. 2000a,b,c; Treu et al. 
2001a, b), are generally carried out under the guiding principle that the FP reflects the 
existence of an underlying mass-luminosity relation for such galaxies (see Faber et al. 
1987; van Albada, Bertin, & Stiavelli 1995), in a scenario where galaxies are homolo- 
gous systems in dynamical equilibrium. According to this simplified picture, galaxies are 
Send offprint requests to: G. Bertin 
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thought to be characterized by a universal internal structure (density and pressure tensor 
distributions), which differs from galaxy to galaxy only in terms of the relevant scales. 
An empirical indication that this may be the case is suggested by the fact that, on the 
large scale, the luminosity profile of elliptical galaxies appears to be universal. 

The empirical i? 1 / 4 luminosity "law" (de Vaucouleurs 1948) has long been recognized 
to fit the surface brightness profiles of elliptical galaxies successfully, to the point that, in 
the absence of other (spectroscopic) indicators, elliptical galaxies are routinely identified 
by means of this characteristic photometric signature. In this sense, it might be argued 
that on the large scale elliptical galaxies are characterized by a universal luminosity profile 
and thus they are also likely to be dynamically homologous systems. In practice, if we 
exclude the central regions of the galaxy, generally affected by seeing or instrumental 
PSF and sometimes studied separately in view of the possible presence of a massive 
black hole, the overall i? 1 / 4 fit is characterized by residuals typically of the order of 
0.1 — 0.2 mag/arc sec 2 (e.g., see de Vaucouleurs & Capaccioli 1979; Capaccioli 1987; de 
Carvalho & da Costa 1988; Capaccioli 1989; Burkert 1993). These deviations from the 
i? 1 / 4 law, although small, are often larger than the typical observational errors involved. 
Interestingly, they are systematic and appear to correlate with the galaxy luminosity 
(Michard 1985; Schombert 1986; see additional comments in the following paragraph). 
In general, the surface brightness profile in the external regions is found to be flatter for 
more luminous galaxies. 

A simple and popular way to model the systematic non-homology noted in the surface 
brightness distribution of early-type galaxies is to consider an i? 1 /™ luminosity profile, 
with n a free parameter, as a generalization of the R 1 ' 4 law (Sersic 1968; see also Davies et 
al. 1988; Capaccioli 1989; Caon, Capaccioli, & D'Onofrio 1993, hereafter CCD93; Young 
& Currie 1994; D'Onofrio, Capaccioli, & Caon 1994, hereafter DCC94; Prugniel & Simien 
1997, hereafter PS97; Wadadekar, Robbason, & Kembhavi 1999). The B}l n "law" has 
been found to be a statistically convenient description to address some issues related 
to the FP of elliptical galaxies (see Graham et al. 1996; Ciotti, Lanzoni, & Renzini 
1996; Graham & Colless 1997; Ciotti & Lanzoni 1997; Graham 1998). In addition, it 
has been applied to the description of the surface brightness profiles of galaxy bulges 
(see Andrcdakis, Peleticr, & Balcells 1995; Courteau, De Jong, & Broeils 1996). On the 
theoretical side, it has been the focus of several general investigations (see, e.g., Makino, 
Akiyama, & Sugimoto 1990; Ciotti 1991; Gerbal et al. 1997; Andredakis 1998; Ciotti & 
Bertin 1999). If applicable, this description would bridge the gap between the classical 
R 1 ' 4 profile and the exponential profile characteristic of disks. It has soon become clear 
that the shape parameter n of the i? 1 /™ law correlates with global quantities such as 
total luminosity and effective radius (CCD93; PS97; see also Khosroshahi et al. 2000). 



4 Bertin et al.: Weak homology of elliptical galaxies 

For example, on a sample of 52 E/SO galaxies with Mb < —17.3, CCD93 find the following 
relation: 

log n = 0.28 + 0.52 log R c (2) 
with rms = 0.18. Here n is the best-fit value obtained by fitting the luminosity profile 
along the major axis, and R e is determined from the curves of growth, independently 
of the fit procedure. CCD93 noted that the transition from normal to bright ellipticals 
occurs at n = 4, R c = 3 kpc, while low-luminosity ellipticals (Davies et al. 1983) are 
characterized by n w 1. The above relation implies the existence of a similar relation 
between n and Mb, since R c and Mb correlate with each other. For example, CCD93 
and DCC94 found: 

n = -19.082 - 1.211Mb ~ -19.082 + 3.0275 logi. (3) 
The average scatter around this relation turns out to be (An) <~ 1.895, but the effective 
scatter increases significantly with L. 

Some confusion has accompanied the application of the R}l n models to observed 
objects. For example, PS97 report the correlations n — R c and n — Mb for a broad 
sample of ellipticals; but here, in contrast with CCD93, the correlations are shown to 
saturate, for logi? [kpc] > 0.8 and Mb < —20, to values of n close to 4. The apparent 
discrepancy is to be traced to the fitting procedure used. In fact, PS97 referred to the 
curves of growth (see comments in Appendix B.l below). As a result, for a given sample 
of early-type galaxies, values of n not larger than 6 have been claimed by PS97, based 
on the use of curves of growth, while values of n up to 15 have been found (see CCD93; 
DCC94) based on direct profile fitting. Curiously, for the low-luminosity objects, where 
the latter procedure sets values of n below 5, the two methods appear to be mutually 
consistent. 

Another natural way to address systematic deviations from the i? 1 / 4 law is to imagine 
(e.g., see Saglia et al. 1997) that the galaxy results from the superposition of an i? 1 / 4 
component (presumably spheroidal) and an exponential component (presumably flat). 
This latter procedure is appealing from the dynamical point of view (SOs are often thought 
to be made like this), yet the procedure is usually applied as a purely photometric method 
(with no check on its dynamical foundation; attempts in this latter direction have been 
made, see Scorza & Bender 1995). 

In view of constructing a physical rationale at the basis of the FP, the existence of 
strict homology is not necessary (see Ciotti & Pellegrini 1992; CCD93; Renzini & Ciotti 
1993; Djorgovski 1995; Capelato, de Carvalho & Carlberg 1995; Hjorth & Madsen 1995; 
Ciotti, Lanzoni, & Renzini 1996; PS97; Pahrc & Djorgovski 1997; and references therein). 
In fact, many of the arguments that allow for a physical interpretation of the observed 
scaling law would carry through even in the case that we may call of weak homology. This 
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is a condition by which the structure and dynamics (density and pressure tensor distri- 
butions) of early-type galaxies vary systematically with galaxy luminosity. An indication 
that this may be a realistic picture is indeed provided by the systematic trends noted 
above in terms of R}l n fits. However, we should emphasize that, based on the photomet- 
ric properties alone of early-type galaxies, we cannot draw conclusions about the FP. For 
this latter goal we have to make use of dynamical models, because we need to translate 
the trends observed in the photometric profiles into trends of the relevant coefficients 
entering the virial theorem. The models used must at least be generally consistent with 
the available kinematical properties of early-type galaxies, otherwise the consequences 
drawn would be based on unjustified models. An interesting attempt at combining pho- 
tometric and kinematic information towards a comprehensive interpretation of the FP 
has been made recently by Gerhard et al. (2001). However, their detailed dynamical 
investigation of a sample of 21 galaxies is far from being conclusive, because their sam- 
ple is mostly biased in the direction of luminous objects and because their kinematical 
measurements, although based on the sophisticated technique of fitting line profiles, are 
admittedly subject to significant model limitations and large uncertainties. 

We should emphasize that observations already provide important empirical evidence 
both against strict homology and against the constant stellar mass-to-light ratio hy- 
pothesis. On the one hand, the existence of a relation such as that given in Eq. (3) and 
the fact that trends in the dark matter content with luminosity are noted readily argue 
against strict homology. On the other hand, the existence of clear correlations involving 
the galaxy luminosity, such as the color-luminosity relation, and the fact that the coef- 
ficients that define the FP depend on the observed waveband suggest that systematic 
changes in the underlying stellar populations are involved. At the interpretation level, 
an explanation based on a pure age-luminosity relation is not viable, because it is hard 
to reconcile with the small and approximately uniform scatter of the FP at zero redshift 
and because it turns out to be incompatible with the observed properties of the FP at 
intermediate redshifts. In addition, it has been shown that interpretations based solely 
on cither a dark matter-luminosity relation or a metallicity-luminosity relation or a sys- 
tematic change of the Initial Mass Function with luminosity are not satisfactory because 
they do not survive a close scrutiny. In this paper we will add further internal consistency 
to the arguments in favor of a form of weak homology at the basis of the FP. 

In Section 2 we carry out an in-depth analysis of 4 significant cases of deviations from 
the R X ' A law by analyzing the surface brightness profiles of NGC 1379, NGC 4374, NGC 
4458, and NGC 4552. We will confirm that a generic i? 1 /™ law, with n a free parameter, 
can provide superior fits (the best-fit value of n can be lower than 2.5 or higher than 
10), better than those that can be obtained by a pure i? 1 / 4 law, by an i? 1 < /4 +exponential 
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model, and by other dynamically justified self-consistent models. Although it may be 
argued that four galaxies are a sample of little statistical significance, the detailed tests 
that we have performed on these otherwise normal objects and, especially, the systematic 
trends with luminosity noted in the literature for large samples (see Eq. (3), based on a 
sample of 52 galaxies) provide convincing evidence that homology cannot be considered 
a strict rule for early-type galaxies. 

In Section 3 we focus on the main implications of these conclusions on the problem of 
the tilt and the thickness of the FP. For a constant (from galaxy to galaxy) stellar mass- 
to-light ratio a suitable correlation n(L) could be consistent with the observed FP. We 
demonstrate that this statement may be true even in the presence of significant scatter 
in the correlation n(L), by means of a general mapping procedure from the space of the 
intrinsic physical model quantities to the space of observed quantities. We also show that, 
in principle, in the absence of strong empirical constraints on the correlation n(L), the 
Fundamental Plane could be realized in an infinite number of ways, two of which are 
indeed the constant stellar mass-to-light solution, and the strict i? 1 / 4 homology solution 
combined with a suitable mass-luminosity relation. Actually, indications provided by 
a set of "data points" associated with a sample of 14 galaxies are shown to suggest 
that neither the strict homology nor the constant stellar mass-to-light solution are a 
satisfactory explanation of the observed Fundamental Plane. 

In the concluding Section 4, we briefly address some interesting consequences on the 
expected trends of the observed FP at intermediate redshifts and outline some natural 
possible extensions of the work presented here. 

Background work is recorded in the Appendices. Appendix A summarizes the prop- 
erties of the models used in the fits presented in Section 2. Appendix B describes the 
relevant fitting procedures. Appendix C illustrates some properties of isotropic dynamical 
models in relation to the definition of the central velocity dispersion cr - Finally, Appendix 
D records the calculation of the virial coefficient associated with the dynamical models 
considered in this paper. 

2. Deviations from the i? 1 / 4 law 

2.1. Preliminaries 

Many studies have already addressed the issue of the universality of the i? 1 ' 4 law (e.g., 
see Burkert 1993; CCD93; DCC94; PS97; and references therein). In this paper we would 
like to focus on some relatively round objects, for which a significant departure from the 
R 1 / 4 law has been noted. These objects will be used as benchmark prototypes on which 
detailed tests will be performed about the quality of a number of photometric models, 
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about various effects associated with the adopted fitting procedure, and about the role 
of the radial range over which the photometric fit is made. 

In general, photometric profiles are compared with the i? 1 / 4 law, or with other pho- 
tometric models, only over a finite radial interval (-R m in, -Rmax)- The inner limit is often 
taken so as to avoid a detailed discussion of the effects of seeing or instrumental PSF (the 
effects of the PSF on the observed surface brightness profile become negligible beyond a 
radial distance from the galaxy center of 2-3 times the scale of the relevant FWHM; e.g., 
see Bertin, Saglia, & Stiavelli 1988); furthermore, the innermost region of a galaxy may 
present features that are related to the presence of a massive central black hole, which 
is generally studied as a separate issue. The outer limit often reflects the quality of the 
available data, since special care has to be taken to properly subtract the contribution 
of the sky surface brightness (special techniques make it possible to measure reliable 
photometric profiles down to [Ib ~ 28 — 29 mag/arc sec 2 ). Of course, a correct estimate 
of the sky background is considerably more difficult in the case of galaxies with nearby 
companions. 

Over a rather small radial range, a very good fit can be obtained rather easily and 
it is thus hard to discriminate among different models. For example, within the radial 
range 0.1i? o < R < l.5R Cl the galaxies NGC 1379 and NGC 4374 are well fitted by the 
i? 1 / 4 law. In fact, based on the definition of average residual 

1/2 

(Aj*> = >i/ 4 (i?i) -Mobs(^)n , (4) 



1 ( N 

* i— 1 



where N is the number of photometric data points fi b s inside the radial range of the fit, 
and of maximum residual 

(A/z) max = max \ni/i{Ri) - Mobs(-R»)| , (5) 

^min <Ri< ^max 

the galaxies NGC 1379 and NGC 4374 are found to be characterized by (A/x) = 
0.033 mag/arcsec 2 , (A/z) max = 0.079 mag/arcsec 2 and by (A/i) = 0.052 mag/arcsec 2 , 
(A^) max = 0.108 mag/arcsec 2 respectively (Burkert 1993). As we will see, these re- 
sults change significantly when the radial range of the fit is increased (but see Makino, 
Akiyama, & Sugimoto 1990). 



2.2. The sample 

In the following we consider four apparently round "prototypical objects", NGC 1379, 
NGC 4374, NGC 4458, and NGC 4552, that are known to present significant deviations 
from the standard R}I A luminosity distribution. We refer to photometric data (in the 
Blue band) available in the literature, taken from Caon, Capaccioli, & Rampazzo (1990) 
and Caon, Capaccioli, & D'Onofrio (1994), that have pointed out the existence of a 
departure from the R 1 ' 4 behavior. Indeed, we have selected the four galaxies, from the 
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Fig. 1. R 1 / 11 fit. R is in arcsec, \i in mag/arc sec 2 . The data points in Figs. 1-5 arc in 
the Blue band, taken from Caon, Capaccioli, & Rampazzo (1990) and Caon, Capaccioli, 
& D'Onofrio (1994). Note the limiting surface brightness reached for each galaxy. 



sample of 80 galaxies studied by these authors, based on the following criteria: (1) Round 
isophotes; (2) Radially extended available photometry; (3) Significant deviations from 
the R 1/4 law. At this stage, we are not interested in the small-scale behavior of the 
luminosity distribution at very small galactocentric radii and so we did not look up 
high resolution data from more modern investigations (e.g., from HST; see Stiavelli et 
al. 2001), which would address issues beyond the goals of this paper. On purpose, the 
adopted criteria l)-3) do not refer to any, possibly available, kinematical information on 
the four objects, because we whish to conform to the usual (i.e., purely photometric) 
selection rules adopted in the construction of early-type galaxy samples for studies of the 
FP. 

For this set of objects we consider the relative performance of four different photo- 
metric models: (i) R}l n (one free dimensionless parameter: n), (ii) i? 1 / 4 (no free param- 
eters), (hi) ii^+exponential (two free parameters: L exp /_L t ot and R eKp / R1/4), and (iv) 
foo models (one free parameter: VP). The main properties of the models are summarized 
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in Appendix A. In all cases, in addition to the best-fit values of the free dimensionless 
parameters, if available, the fitting procedure leads to the determination of two scales, 
i.e., of the model effective radius R c and of the model total luminosity L. 

Each of the four different photometric models is tested against a variety of options for 
the fitting procedure: (a) profile fitting or circularized profile fitting or curve of growth 
fitting, (b) reduced or extended radial range, and (c) constrained or unconstrained inte- 
grated luminosity. These fitting procedures are briefly described in Appendix B. 

2. 3. Relative performance of four photometric models for four prototypical objects 




1 1.5 2 2.5 0.5 1 1.5 2 2.5 

log(R) log(R) 

Fig. 2. R 1 / 11 fit. R is in arcsec, fi in mag/arc sec 2 . 



In Table 1 we record the results for the best-fit i? 1 /™ models, for the reference option 
of using circularized profiles, of considering the most extended radial range available, 
and of imposing no constraint on the integrated luminosity. The quality of the fit is 
measured by the parameter % 2 , defined in Appendix B, and is also described here by the 
values of (A^t) and (A/i) max (sec Eqs. [4]-[5]). Some of the digits of the various entries 
are not significant, but we have decided to keep them since there should be no reason for 
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misinterpretation. The last two columns record the surface brightness excursion available 
for the photometric data and the total absolute magnitude of each object. The best-fit 
R l / n models are then illustrated in Figs. 1 and 2. 

Similarly, in Table 2 we summarize the corresponding performance of the other three 
photometric models (for the same option of fitting procedure as in Table 1). It is clear 
that, especially for the two objects (NGC 4374 and NGC 4552) that are characterized by 
high values of n, the model effective radii thus determined by the fit can vary dramatically 
(but see comments in Sect. 2.2). The fits are then illustrated in Figs. 3-5. For the low- 
n galaxy NGC 1379, the models are not providing a satisfactory fit, which is not 
surprising since for this galaxy the role of rotation is likely to be important (see the 
kinematical profile recorded by Graham et al. 1998) and thus the galaxy falls outside the 
physical scenario at the basis of the models. In addition, two minima in x 2 arc actually 
available for the foo models for this galaxy: the higher \ 2 model (lower not illustrated 
in Fig. 5) is less bright than the galaxy in the central region (opposite to the behavior 
of the model fit shown in the Figure), both models are brighter than the galaxy at large 
radii. 

In Table 3 we summarize the results of the best-fit i? 1 /™ models obtained under 
various fitting options. It is clear that many of the quantities thus derived are sensitive 
to the adopted fitting procedure. As an extreme example, we point to the case of NGC 
4374, for which the best-fit model for the reduced radial range (case b of Table 3) would 
be characterized by R c — 1092 arcsec, so that most of the total luminosity of the galaxy 
would be assigned to the unobserved part of the photometric profile. Our results quantify 
and generally confirm trends already noted in earlier studies (see Introduction; but see 
also some important remarks in Appendix B.l). 

Table 1. R}l n best fits. Circularized profiles, extended radial range, no L constraint. R c 
is expressed in arcsec. 



NGC 


x 2 


Re 


0FP 


n 


(A M > 


(A/i) max 


/J-max /J-min 


Ml 


1379 


0.31 


25.7 


-5.33 


2.38 


0.042 


0.086 


7.64 


-19.51 


4374 


0.46 


259.0 


-5.08 


9.48 


0.044 


0.123 


7.54 


-21.35 


4458 


0.85 


22.2 


-5.57 


3.19 


0.044 


0.115 


4.76 


-18.45 


4552 


1.32 


104.0 


-5.13 


11.14 


0.095 


0.239 


8.30 


-20.73 



a 0pp is given by Eq. (6), with R e in arcsec and SB e in mag/arc sec 2 . 

b absolute magnitude in the Blue band (apparent magnitude and extinction from de Vau- 
couleurs et al. (1991), adopted distance of 18.3 Mpc) 
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Table 2. Comparison of best-fit models. Circularized profiles, extended radial range, no 
L constraint. Columns from second to fourth refer to i? 1 / 4 models, from fifth to ninth to 
i? 1 / 4 +exponential, and last four to f x models. R e is expressed in arcsec. 



NGC 


x 2 


Rc 


0FP 


x 2 


Rc 


0FP 


Lcxp/ Ltot 


iW-Ri/4 


x 2 


R c 


4>FP 




1379 


5.85 


34.9 


-5.33 


0.78 


27.9 


-5.31 


0.22 


0.65 


3.22 


32.6 


-5.30 


9.6 




















(5.22) 


(46.9) 


(-5.38) 


(5.2) 


4374 


9.05 


75.2 


-4.96 


0.80 


104.0 


-4.99 


0.06 


0.06 


1.56 


98.5 


-4.94 


7.3 


4458 


1.13 


26.6 


-5.57 


0.74 


21.4 


-5.55 


0.13 


0.96 


0.71 


23.9 


-5.53 


9 


4552 


12.23 


41.2 


-5.07 


3.07 


58.1 


-5.10 


0.10 


0.09 


2.86 


146.4 


-5.22 


6.3 



a effective radius of the combined profile 



Table 3. Best-fit R}l n models under various fitting options. R c is expressed in arcsec. 



NGC 


n a 


Rc 


0PP 


x 2 


n b 


Rc 


0FP 


x 2 


n c 


Rc 


0FP 


1379 


2.88 


26.8 


-5.32 


0.38 


2.36 


25.4 


-5.33 


0.31 


2.39 


25.7 


-5.33 


4374 


5.73 


124.1 


-5.00 


0.33 


14.5 d 


1092 


-5.25 


0.45 


9.31 


247.0 


-5.07 


4458 


4.61 


29.8 


-5.57 


0.83 


3.28 


22.7 


-5.57 


0.81 


3.21 


22.3 


-5.57 


4552 


5.87 


56.7 


-5.08 


0.85 


9.95 


95.2 


-5.12 


1.33 


11.9 


119.9 


-5.15 



a extended radial range, curves of growth, no L constraint (constant photometric error as- 
sumed, see App. B.l) 

b reduced radial range (/i ma x — /imin =4.5 mag), circularized profile, no L constraint 
c extended radial range, constrained luminosity L, circularized profile 
d solution "crossing" (see description in App. B.3) 



For the set of four galaxies, considered to be as prototypes of objects with a genuine 
deviation from the R 1 ' 4 law, we have found the following results. These are based on 
direct inspection of Figs. 1-5 and especially on comparison of the x 2 values relative to 
the various fitting options in Table 1 and 2. For all galaxies the fit by a pure i? 1 / 4 law is 
the worst, with significant residuals. For one galaxy (NGC 4458) the i? 1 / 4 +exponential 
model gives a marginally better fit with respect to the R 1 / 71 model (x 2 = 0.74 vs. x 2 = 
0.85), but for this galaxy also the pure i? 1 / 4 law and the foo models give reasonable 
fits (x 2 = 1-13 and 0.71, respectively); not surprisingly, this is the galaxy with the 
smallest surface brightness excursion (^ max — /i m in) sampled by the available data (less 
than 5 magnitudes) and the one for which the best-fit value of n (n w 3.2) is closest 
to the standard value of 4. For the three remaining galaxies (NGC 1379, NGC 4374, 
NGC 4552), for which the surface brightness profiles extend over a range of more than 7 
magnitudes, the R}l n model definitely gives the best fit, even if it has fewer parameters 
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than the i? 1 / 4 +exponential model. In addition, the two objects (NGC 4374 and NGC 
4552) that turn out to be characterized by high values of n (of order 10), when fitted 
by the i? 1 / 4 +exponential model, would lead to a disk component with an unusually 
high central brightness (18.2 and 18.1 mag/arcsec 2 in the Blue band, respectively), at 
variance with the central brightness of galaxy disks that characterize spiral galaxies. It 
would be interesting to check wheter this is a common feature of i? 1 / 4 + exponential 
decompositions of early-type galaxies. 

2.4- Robustness of a photometric indicator 

In the discussion of the properties of the FP, it has often been noted (see Kelson et al. 
2000a; Treu et al. 2001a; see also Trujillo, Graham, & Caon 2001) that much of our 
uncertainties in relation to the parameters derived by fitting the photometric profile is 
resolved by an unexpected circumstance, that is the empirical robustness of the photo- 
metric indicator that enters the FP with respect to the modeling procedure. The relevant 
indicator is 

(hp = log R C -/3SB C , (6) 

with (3 « 0.32 (see Eq. (1) and following remarks). 

An important point should be noted here. There are two different ways of considering 
possible variations of </>fp for a given object, on the basis of the definition of SB C reported 
just before Eq. (1). One is that of considering both R e and L as model quantities, both 
determined by the fit to the data (see the comments on the extreme case of NGC 4374 
made at the end of the previous subsection). The other is to keep the total luminosity 
as fixed (for the set of four galaxies studied earlier in this paper, see the last column of 
Table 1), and let only R e be determined by the fit. In this latter case, it is readily shown 
that the change in the value of the photometric indicator associated with a small change 
of R c due to the modeling procedure is <50fp = (1 — 5/3) log e 5R e /R e ~ —0.26 SR e /R e ; 
this shows that the change in (/>pp is, in this case, systematically smaller than the relative 
change in R c . No matter which definition of 54>fp we refer to, the impact of a variation 
in (j)pp on the FP itself can be quantified by stating that, if we call 1Z the ratio between 
the value of R e predicted by the FP, at a given value of do, and the observed value, then 
5K/K = -(50 FP /loge w -2.3 <ty FP . 

Here, based on the first definition of SB e (for which L is taken as a model quantity), 
we confirm the general validity of the statement. Surprisingly, the photometric indicator 
is robust both with respect to changes of the adopted photometric model (compare 
Table 1 and Table 2) and with respect to variations in the fitting procedure within the 
same photometric model (compare Table 1 and Table 3). Therefore, for the purpose of 
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measuring the quantities entering the FP we are on safe grounds, whereas for the purpose 
of interpreting the FP we are left in a rather ambiguous situation. 
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2.5. The kinematical counterpart. Low-luminosity galaxies as a physically separate class 
of objects 

The investigation provided earlier in this paper sheds some light on the issues of modeling 
and confirms that homology cannot be considered a strict rule for early-type galaxies. 
However, the discussion so far has been restricted to the photometric aspects, while 
the possibility of weak homology (i.e., of systematic trends in structural and dynamical 
properties as a function, e.g., of galaxy luminosity) should be assessed by considering the 
entire set of quantities that enter the FP. In particular, it is unavoidable that one should 
address the internal kinematical properties of the galaxies involved. In fact, it has long 
been known that low-luminosity ellipticals are likely to be rotation supported systems 
(Davies et al. 1983). If we are interested in drawing a connection between the low values 
of n that best describe low-luminosity ellipticals and the dynamical coefficient that enters 
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the virial theorem, we should not forget that the virial theorem relation reflects the entire 
dynamical structure of the system involved. 

To be more specific, it is often argued that the relevant dynamical coefficient can 
be conveniently calculated under the simplifying assumption that the system be char- 
acterized by an i? 1 /" profile in projection, while intrinsically the system be spherical, 
non-rotating, and strictly isotropic. From the Jeans equations such coefficient is deter- 
mined in a straightforward manner (see PS97; Ciotti & Lanzoni 1997; see also Sect. 3 
and Appendix D below). In practice, it remains quite doubtful that such coefficient can 
be usefully inserted into a realistic discussion of the FP, given the fact that the velocity 
dispersion profiles thus calculated are often qualitatively different from those observed 
(see Appendix C). Indeed, some of the difficulties that we may find in placing the four 
galaxies described in this section into the context of the FP, as will be attempted in 
Sect. 3 (in particular, see Figs. 6-9), may just be due to this basic modeling problem. 
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So far we may claim that we have a reasonable representation of the dynamical 
structure only of bright ellipticals, at least in the case of intrinsically fairly round objects 
(e.g., see Bertin & Stiavelli 1993). Low-luminosity ellipticals appear to be a physically 
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separate class of objects. Until we gain, by systematic spectro-photometric surveys and 
by new and accurate modeling, a good picture of how early-type galaxies are structured, 
from faint to bright objects, there is no hope that we can make a reliable discussion 
of the impact of the observed deviations from the i? 1 / 4 law on the FP. Nonetheless, in 
the absence of such comprehensive dynamical modeling, by using the one-component, 
spherical, non-rotating, isotropic i? 1 /™ models (even if wc know them to be an inadequate 
representation) we will be able to demonstrate, in the next Section, that the absence of 
strict i? 1 / 4 homology can be consistent with an entire spectrum of possible interpretations 
of the FP. Curiously, these interpretations may well hold even if the empirical relation 
n = n(L) is affected by large scatter. 



3. Impact on the Fundamental Plane 

In order to assess the consequences of weak homology on the interpretation of the FP, 
below wc follow a line of thought already introduced, e.g., by Ciotti & Pellegrini (1992). 

For a stationary stellar system, the scalar virial theorem can be written as 
GT~L 



n = Kv< . (?) 

He 

where the coefficient Ky takes into account projection effects, the specific star orbit 
distribution and the specific density distribution, and effects related to the presence 
of dark matter; here Y* is the stellar mass-to-light ratio in the given band used for 
the determination of L and R c . Equations (1) and (7) imply that in real galaxies, no 
matter how complex their structure is, the dimensionless factor Y*/Ky is a well-defined 
function of any two of the three observables [L, R e , a®). For example, by eliminating <j$ 
from Eqs. (1) and (7) we obtain 

— OC H ^-WP+a)/a L (5P-a)/ a _ ^ 

K\ 

This fact is remarkable by itself. All stellar systems described by Eq. (7) are in virial 
equilibrium, but only those stellar systems for which Y*/Ky scales according to Eq. (8) 
are placed onto the observational manifold given by Eq. (1). Of course, this requirement 
is only necessary for a model to describe real galaxies. In fact, the FP is not uniformly 
populated by galaxies; instead, observed stellar systems populate only a well-defined 
region of the FP. 

After recalling two simple extreme interpretations of the FP (Subsection 3.1), in 
Subsection 3.2 we move to a more general exploration, which will bring out once more a 
problem of fine tuning met in the study of galaxies. 
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3.1. Two extreme interpretations of the Fundamental Plane 

The possibility of simple extreme explanations of the FP is associated with the fact that 
the exponent of the effective radius appearing in Eq. (8) is observationally very small, 
i.e., (2 - 10/3 + a) I a w 0.04 in the Blue band. 

A first possibility to reproduce the scaling relation in Eq. (8) is then to assume 
homology (i.e., Ky identical for all galaxies), a variable stellar mass-to-light ratio (see, 
e.g., Renzini & Ciotti 1993; van Albada, Bertin, & Stiavelli 1995), and the empirically 
suggested identity 2 — 10/3 + a = 0. With these assumptions we find 
T, oc L* ; (9) 
here 6 = (2 — a) /2a ~ 0.30 ± 0.064, where the scatter has been calculated by standard 
error propagation techniques. From a physical point of view, in a simple (coeval and 
single metallicity) monolithic formation, passive evolution scenario, a mass-luminosity 
relation for galaxies, of the kind suggested by Eq. (9), can result if, with significant fine 
tuning, the basic properties of the Initial Mass Function, which describes the origin of 
stellar populations, properly correlate with galaxy mass. 
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An alternative extreme explanation can be proposed by assuming a constant stellar 
mass-to-light ratio and the existence of weak homology (see Ciotti & Pellegrini 1992; 
CCD93; Rcnzini & Ciotti 1993; Djorgovski 1995; Hjorth & Madsen 1995; Ciotti, Lanzoni, 
& Renzini 1996; PS97; and references therein). In this case it is the quantity Ky (instead 
of T*) that is required to be a well-defined function of R c and L. If we take 2— 10/3+a = 0, 
the required dependence is 

K v oc L~ s , (10) 
with S and its scatter the same as above. As noted in the Introduction, this latter interpre- 
tation is qualitatively encouraged by the observations. Quantitatively, we can easily show 
that we are in the right ballpark. In fact, from Eq. (10) it is realized that a factor of 20 
in L, would require a change in Ky by a factor of sa 2.45. If we consider one-component, 
spherical, non-rotating, isotropic R}l n models and a simulated spectroscopic aperture of 
Rc/8, the relation Ky(n) can be computed (Ciotti & Lanzoni 1997; PS97; see Appendix 
D). An accurate and convenient analytical representation in the range 1 < n < 10 (with 

typical errors on the order of one percent) is given by 
73 32 

K ^ - 10.465 + (n-0.94) 2 +0 - 954 - ^ 
It is then demonstrated that the observed excursion in n (see Sect. 2) could well be 

consistent with the desired range in Ky. A similar argument also holds for the 

models, for which we have the following simple analytical interpolation formula, accurate 

to around one percent, 

T . , T . 142.3- 41.51* + 2.66* 2 

Ky(^) ~ = (12) 

v ; 30.61 - 10.7* + * 2 
in the range 2 < * < 10. 

From this simple discussion we might conclude that weak homology of elliptical galax- 
ies might indeed be consistent with a constant T* and the FP. 

Unfortunately, the determination of the coefficient Ky is sensitive to the choice of 
models (see the expressions recorded above and App. D) and to the choice of the aper- 
ture radius. For the R}' n models, the latter dependence is not too significant, even in 
the presence of the well-known unrealistic central depression in their projected velocity 
dispersion profiles (see App. C); this has been checked by comparing the values of Ky 
obtained on the basis of an i? e /10 aperture with those based on Re/8. In order to make 
a physically better justified case, one might think it more appropriate to refer to the /oo 
models, which are known to possess realistic kinematical profiles. Note that this would 
lead to significantly different answers with respect to the application of R}' n models. For 
example (see App. D), if we take NGC 4552, the i? 1 /™ modeling would give n ~ 11.14 and 
Ky ~ 1.6, while the /oo modeling would set * ~ 6.2 and Ky ~ 2.5; for the galaxy NGC 
4458 (closest to the standard i? 1 / 4 law in our sample), we find n ~ 3.2 and Ky ~ 5.7 
based on i? 1 /™ modeling or * ~ 9, Ky ~ 3 based on modeling. Here we prefer not 
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to pursue this line of investigation, because the empirical relation = ^(L) has not yet 
been studied systematically and because, in any case, we know a priori that such model- 
ing, even if it might turn out to be partially successful, would lack physical justification 
when applied to low-luminosity, rotation-supported galaxies. 

3.2. A more general interpretation allowed by weak homology 

We now describe a more general interpretation of the FP allowed by the existence of 
weak homology. For simplicity, we refer here to the family of R}l n models that has led to 
Eq. (11), although the method proposed could be implemented for any well-defined family 
of dynamical models. In particular, here we explore the theoretical manifold obtained by 
mapping, with the aid of Monte Carlo simulations, the model space (n, T_, L, R c ) into the 
observed space (L, i? , (T ). The domain of model galaxies considered in the simulations is 
defined by 1 < R c < 20, 0.1 < L < 50, and 100 < Co ^ 350, where radii are measured in 
kpc, luminosities in 10 10 solar Blue luminosities, and velocities in km s -1 . For each four- 
dimensional point (n, T*, L, R e ) in model space we calculate cr from Eq. (7) and then 
we check if the mapped (virialized) system belongs to the FP. This allows us to select in 
model space, as acceptable candidates for real galaxies, only those points that turn out 
to be compatible, in the observed space, with the FP correlation and the observed scatter 
around Eq. (1). Wc also explore the "backbone" model manifold obtained by artificially 
reducing the FP thickness to a limiting zero-scatter case. For simplicity, we plot the 
projections of the manifold thus derived on the natural planes (L, Y*) and (L, n), which 
are precisely the two planes that are used to illustrate the two extreme interpretations 
described in Subsection 3.1. 

We start this exploration by applying the proposed technique to the two cases dis- 
cussed in the previous subsection. The case of R 1 / 4 homology is studied by restricting the 
model population to a set of points for which 3.99 < n < 4.01. The derived dependence 
of Y» on L is well fitted by the power law relation of Eq. (9), with 5 ~ 0.27 (consistent 
with the value 0.30± 0.064 argued in Subsection 3.1), with reduced x 2 — 1-7. The scatter 
of points along the vertical axis at fixed L (see Fig. 6) reflects the assumed thickness of 
the FP. Indeed, if we artificially reduce the scatter in R e from 0.15 to 0.01, the derived 
power-law exponent becomes 5 = 0.28, with reduced x 2 ~ 7 x 10~ 2 (the thick locus of 
points in Fig. 6). 

To illustrate the second case (see Fig. 7), we pick a constant stellar mass-to-light ratio 
3.49 < Y» < 3.51. The selected points are quite dispersed along the n axis. A good fit is 
obtained by a law of the form n = a + 61og_L, with b ~ 4.02 and \ 2 ~ 1.3, not so well in 
agreement with the observed scaling relation given in Eq. (3). After artificially reducing 
the scatter around the FP to 1 percent in i? , we now have b ~ 5.03 and x 2 ~ 0.24. 
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Fig. 6. Projection of the manifold in model space identified by i? 1 / 4 homology. The 
scatter of points reflects the adopted scatter around the FP in the observed space. The 
symbols that trace the thick line represent the "backbone" manifold. Open squares rep- 



and Dl; crosses represent the same four galaxies, following the modeling summarized 
in Tables 2 and Dl (see Subsect. 3.3 for additional description). Filled squares represent 
the sample of ten galaxies considered in Table 5 of Bertin ct al. (1994); see Subsect. 3.4 
of the present paper for a description. 

(By increasing or decreasing the reference value of T*, the simulation points would move 
below or above the strip identified by the present choice of T*.) 

A somewhat unexpected result of these simple simulations is that the relatively small 
dispersion of R c for real galaxies, around the value associated with the FP at any fixed L 
and (To , is consistent with a significantly larger dispersion of the model parameters along 
the vertical axes in Figs. (6) and (7) at fixed values of L. Therefore, weak homology (e.g., 
a correlation between the structural parameter n and the luminosity L) can be at the 
origin of a relatively thin FP, even in the presence of a substantial dispersion of n at 
given L. 

At this point we can proceed to consider the most general case, where neither ho- 
mology is strict nor mass-to-light ratio is constant. In fact, we leave the model points 
to run over the entire ranges of T* values and n values in the domain indicated at the 
beginning of this subsection. For simplicity, we look for the backbone manifold selected 
by artificially reducing the scatter on the observed R e to 0.01. In this way the scatter 
presented by the points on the two natural planes (L, Y») and (L, n) will be determined 
only by projection effects (of the relevant manifold onto the desired planes). 
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Fig. 7. Projection on the plane (L,n) of the manifold in model space identified by the 
constant mass-to-light ratio assumption ~ 3.5. The scatter of points reflects the 
adopted scatter around the FP in the observed space. Solid symbols represent the "back- 
bone" manifold. Open squares represent the four galaxies of Sect. 2, following the i? 1 /™ 
modeling summarized in Tables f and Df (see Subsect. 3.3 for additional description). 

In Fig. 8 we show the projection of what we have called the backbone manifold onto 
the (L, T*) plane. Note how, at any given value of L, a larger range of Y* values is 
permitted, even if the FP considered in the simulations is extremely thin. Therefore, in 
the general case, T* is not necessarily well correlated with galaxy luminosity, as a result of 
the presence of the free parameter n. It is interesting to check how models with different n 
are distributed inside this projection of the backbone manifold; the open circles represent 
models characterized by 3.8 < n < 4.8, a range slightly broader than that shown in the 
homology case shown in Fig. 6. In any case, also under the more general circumstances 
considered here, the FP relation is associated with an overall increasing trend for the 
stellar mass-to-light ratio T* with L. 

Similar comments apply to the projection illustrated in Fig. 9. Here the open circles 
represent models with the stellar mass-to-light ratio T* w 3.5. Note the large spread in 
n present, still consistent with the existence of a thin FP, once the restriction of nearly 
constant T* is relaxed. Note also that there remains only a mild trend of n with L 
required by the existence of the FP, but only if we focus on the zero-scatter limit. Since 
in reality some correlation is observed (see Introduction and Eq. (3)), we should conclude 
that something else, in addition to the FP, characterizes the underlying weak homology 
of real galaxies. 
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Fig. 8. Projection on the (L, Y*) plane of the entire backbone manifold identified by 
the FP with scatter artificially reduced to 0.01. Open circles correspond to models with 
3-8 < n < 4.2. Open squares, crosses, and filled squares define observed galaxies as in 
Fig. 6. 
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Fig. 9. Projection on the (L,n) plane of the entire backbone manifold identified by the 
FP with scatter artificially reduced to 0.01. Open circles correspond to models with 
3.4 < T* < 3.6. Open squares define observed galaxies as in Fig. 7. 



3.3. Locating the four prototypical objects without detailed stellar dynamical modeling 

With the aim of a better appreciation of the naive indications provided by a pure photo- 
metric modeling we have recorded, in Figs. 6-9, as open squares (labels a, b, c, d following 
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the order of Table 1), the points available for the set of four galaxies discussed in the 
previous part of the paper and from the available data about their central velocity dis- 
persions. This information (see Table D.l), based on the i? 1 /™ modeling, allows us to 
extract values of from the observed values of a (119, 296, 106, and 269 km s _1 for 
NGC 1379, 4374, 4458, and 4552, respectively, McElroy 1995). In Fig. 6 and Fig. 8 we 
also record the different points that would be argued on the basis of a one- component 
/oo photometric modeling. 

For the purposes of this Section, the relevant quantities associated with the four 
galaxies are rescaled, when necessary, to the long cosmological distance scale. This is 
consistent with the adopted choice of coefficients for the FP (see Eq. (1) and following 
comments) and with the analysis to be presented in the following Subsection. 

3.4. Trends from a sample of relatively bright galaxies and the role of dark matter 

In order to better relate the arguments provided in this Section with the trends observed 
in real galaxies, we have included in Figs. 6 and 8 some "data-points" that are external 
to the simulations described above. These data-points have been obtained as follows. 

We have referred back to the results of an ESO Key Programme "A search for dark 
matter in elliptical galaxies" and, in particular, to a set of 10 galaxies the properties of 
which are summarized in Table 5 of Bertin et al. (1994). The objects are: NGC 1399, 
NGC 3379, NGC 4374, NGC 4472, NGC 5812, NGC 5813, NGC 7507, NGC 7626, NGC 
7796, and IC 4296 (some of these are also part of the sample considered by Gerhard et 
al. 2001). The numbers from 1 to 10 of the filled squares in Fig. 6 and Fig. 8 represent 
precisely these objects in the order followed above. These galaxies had been considered 
as those for which extensive and reliable stellar dynamical modeling was available, on the 
basis of accurate and radially extended stellar dynamical data. We borrow the values of 

and R e listed in Table 5 of Bertin et al. (1994) under the columns (Ml/Lb)2C and 
R™, respectively, and thus obtain the values of Ky for those objects, based on a detailed 
stellar dynamical modeling. One object, NGC 4374, is also part of the set of four galaxies 
for which we have re-examined the photometric profile in Sect. 2. 

We note that those models did include significant amounts of dark matter. Therefore, 
it is not surprising to find now that several filled squares fall outside the strip outlined 
in Fig. 6, which has been calculated on the basis of strict i? 1 / 4 homology. Then Fig. 8 
shows that the same data-points are well consistent with a possible weak homology that 
is neither the strict i? 1 / 4 homology nor the one based on a constant T*. Note that the 
broad area outlined in Fig. 8 has been obtained by reducing artificially the scatter of the 
FP to 0.01. 
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4. Conclusions 

This paper, organized in two parts, has led to a number of results, although many are 
the questions that remain only partly answered. Here we provide a summary of the main 
points made in the paper and of the most important questions that are left for future 
investigations. 

In the first part of the paper we have re-examined the surface brightness profiles of 
NGC 1379, NGC 4374, NGC 4458, and NGC 4552, four galaxies that are known to present 
significant deviations from the i? 1 / 4 law. We have considered the relative performance 
of four different photometric models, namely the i? 1 /™, i? 1 / 4 , i? 1 / 4 +exponential, and f x 
models. Each of the four different photometric models has been tested against a variety 
of options for the fitting procedure. We have thus reached the conclusion that there is 
indeed convincing evidence that homology cannot be considered a strict rule for early 
type galaxies. In particular, we have found that: 

— For the four galaxies of this limited sample the fit by a pure i? 1 / 4 law is the worst. 

— For three objects the R}l n photometric model definitely provides the best fit, even 
better than that by the i? 1 / 4 +exponential model, although the latter model has the 
advantage of one additional free parameter. 

— For NGC 4458, which is the object with the smallest surface brightness excursion 
available (less than five magnitudes), the i? 1 / 4 +exponential model gives a marginally 
better fit with respect to the i? 1 /™ model; but, for this, the i? 1 / 4 and the f x models 
give practically equivalent fits. 

— The /oo models manage to provide a reasonable description of the bright ellipticals 
even when the i? 1 /™ modeling suggests a value of n significantly larger than 4. How- 
ever, they do not appear to be viable for the low-luminosity ellipticals, which arc 
characterized by a small value of n. This does not come as a surprise, because the 
physical arguments that have led to the construction of the /oo models reflect the pic- 
ture of low rotation, collisionlcss collapse, while low-luminosity ellipticals are known 
to be generally rotation supported. 

— The values of the effective radius resulting from the fit can change dramatically (by 
factors of up to three) , depending on the fitting procedure and on the adopted pho- 
tometric model. On the other hand, we confirm that the determination of the pho- 
tometric indicator, i.e. the combination of photometric quantities entering the FP, is 
robust. 

The good performance of the R}/ n model in describing the photometric profiles of 
galaxies deviating from the i? 1 / 4 law and the existence of an empirical correlation between 
n and L (see Eq. (3), based on a sample of 52 galaxies) suggest that elliptical galaxies 
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might be considered as weakly homologous systems. Of course, a clear indication of 
weak homology cannot be based on photometric properties alone. In the second part 
of the paper, we have then addressed the impact of weak homology on the physical 
interpretation of the FP, starting from a re-discussion of the rationale behind two extreme 
interpretations of the tilt of the FP. One is obtained by assuming strict homology and a 
suitable relation between the stellar mass-to-light ratio and the galaxy luminosity. The 
other is obtained by taking a constant stellar mass-to-light ratio (from galaxy to galaxy) 
and a suitable form of weak homology. The latter interpretation usually raises concerns 
(besides the problems related to the interpretation of the observed correlations involving 
luminosity, such as the color-luminosity relation; see discussion in the Introduction), 
because the large scatter around the correlation n = n{L) appears to be hard to reconcile 
with the small dispersion observed around the FP. We have thus obtained the following 
results: 



— In order to address the problem in its generality, we have developed a systematic 
approach to the study of the FP, based on a mapping between model and observed 
quantities. For simplicity, we have used, as a guiding tool for the present investiga- 
tion, simple one-component, spherically symmetric, isotropic R}l n models. The study 
would be best carried out with other families of models, better justified from the phys- 
ical point of view, but for the present purposes the simple B}l n models used provide 
an adequate demonstration. The use of one-component models is equivalent to assum- 
ing that dark matter, if present, is distributed in direct proportion to the luminous 
matter. 

— Within this framework, we have shown how the scatter of the stellar mass-to-light 
ratio (in the first of the two simplest interpretations) or of the index n (in the second 
interpretation) is fully controlled by the observed thickness of the FP. This quantifies 
in detail the problem of parameter fine tuning that has to be faced. 

— We have then relaxed the hypothesis of a simple (i.e., one-parameter) origin of the 
tilt of the FP and focused, for simplicity, on the structure of the FP with artificially 
reduced thickness. A fine tuning problem obviously remains, because the existence of 
the FP imposes a strict condition on Y*/Ky. However, we have been able to demon- 
strate that the ranges of variation for the stellar mass to light ratio and the index n 
of the models, at fixed galaxy luminosity, can be substantial and still be consistent 
with an arbitrarily thin FP. The fine tuning involved requires that variations in the 
stellar mass-to-light ratio be compensated for by corresponding variations in galaxy 
structure. 
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— Preliminary indications provided by a set of "data points" associated with a sample 
of 14 galaxies suggest that neither the strict homology nor the constant stellar mass- 
to-light solution are a satisfactory explanation of the observed Fundamental Plane. 

This paper opens the way to the following new investigations, to be carried out in 
separate papers: 

— A study based on the approach introduced in Sect. 3 within a more complex dynamical 
structure (presence of dark matter and more complicated geometries for the model 
galaxies) . 

— A study of the evolution of the FP with rcdshift, within a physical interpretation 
based on weak homology. If the observations were able to tell the time evolution of 
the coefficients a and in the exponents in Eq. (8), then, from our mapping proce- 
dure, we could infer the actual evolution shared between luminosity and structural 
parameters. From a different, predictive point of view, if we argue that the tilt of the 
FP is totally due to weak homology (with constant T* from galaxy to galaxy, so that 
the FP is generated by systematic structural variations with luminosity of the quan- 
tity Ky), then we would expect only a parallel shift of the FP, i.e. only a systematic 
shift of its zero-point 7 with time. Clearly, this picture must be only approximately 
viable, because galaxies follow the metallicity - velocity dispersion relation (see Ben- 
der, Burstein, & Faber 1993 and references therein) and the color-magnitude relation 
(see Bower, Lucey, & Ellis 1992), consistent with the fact that the tilt of the FP 
depends on waveband (e.g., see the results based on near-infrared data obtained by 
Pahre, Djorgovski, & de Carvalho 1998). On the other hand, if, as is frequently done, 
we refer to the picture where all galaxies are strictly homologous systems, so that the 
tilt of the FP is completely due to a luminosity dependence of T* (at fixed metal- 
licity), then the luminosity evolution of faint galaxies should be faster than that of 
bright galaxies, and the tilt of the FP would be expected to change with z, with the 
larger changes at the FP faint end. Here several different physical scenarios would be 
available (see also Renzini & Ciotti 1993). 

— A study of the relation between weak homology and the popular paradigm of galaxy 
merging. 

From the observational point of view, it is especially urgent to carry out a discussion 
of weak homology, extended to the spectroscopic (stellar-dynamical) data, for a complete 
sample of nearby early- type galaxies. 
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Appendix A: Photometric Models 

A.l. B}l n and R 1 / 4 models 

Consider the surface brightness profile (Sersic 1968) given by 

I(R) = Joe-"" 17 " , (A.l) 
where 77 = R/R c , n is a positive real number, and b — b(n) a dimensionlcss constant 
such that i? e is the effective radius, i.e., the projected radius encircling half of the total 
luminosity associated with I(R). As proved in a separate paper (Ciotti & Bertin 1999), 
the first terms of the asymptotic expansion of b(n), for large n, are b(n) ~ 2n — 1/3 + 
4/(405n), already accurate to a part in one thousand for n as low as unity. The proper 
normalization of the R 1 /™ profile, to be considered for a case with given scales L and R e , 
is 

where V represents the complete Gamma function. The photometric profile is thus charac- 
terized by two scales (e.g., L and R c ) and one dimensionlcss parameter (n). The standard 
de Vaucouleurs (1948) profile is obtained by setting n = 4. 

A. 2. Two-component R X I A +exponential models 

For the surface brightness profile we take here the superposition of an i? 1 / 4 component 
and a standard "exponential disk" , i.e. 

I(R) = J 1/4 (0) cxp[-6(4)(i?/i? 1/4 ) 1 / 4 ] + 7 exp (0) cxp[-6(l)(i?/i? cxp )], (A.3) 
where R1/4, and i? cxp are the effective radii of the i? 1 / 4 and of the exponential component, 
respectively. By fitting this model to an observed photometric profile, we can derive the 
total luminosities L C x P , £1/4 associated with the separate components and from these 
the total luminosity L tot = Z/1/4 + i cxp - In conclusion, this two-component photometric 
model is characterized by two scales (Ltot and R c ; for the latter quantity we refer to 
the effective radius for the combined profile) and by two dimensionless parameters (e.g., 
Lexp/Ltot and Re^-p / R\ / 4) ■ 

A.3. foo models 

The /oo models are dynamical models applicable to collisionless stellar systems. They 

have been constructed (Bertin & Stiavelli 1984) by following the physical scenario that 

elliptical galaxies may have formed through collisionless collapse (van Albada 1982). In 

the spherical limit, their anisotropic distribution function is given by 
J A(-Er/ 2 cxp(-aE-cJ 2 /2) if E < 0, 

J co — \ (A-4) 

{ i£E>0, 
where E — v 2 /2 + <&(r) and J are the star energy and the magnitude of the star an- 
gular momentum, per unit mass, respectively. Two of the three available dimensional 
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constants (A, a, and c) set the physical scales of the model (e.g., the total mass and the 
half-mass radius Tm)- The remaining dimensionless parameter, for example written as 
7 = ac/(4irGA), then characterizes the one-parameter family of foo models. The solution 
of the Poisson equation provides the relation 7 = 7(*), where denotes the dimension- 
less central potential W — — a$(0). In practice, under the assumption that the stellar 
mass-to-light ratio is constant within the galaxy and after projection for comparison 
with the observations, the relevant photometric profile associated with the models, 
in addition to a luminosity scale (L) and a length scale (R e w rAf/1.3), depends on one 
dimensionless parameter only (\D r ). A detailed discussion of the properties of this equilib- 
rium sequence is given by Bertin & Stiavelli (1993). Here we limit our discussion to 
models characterized by 2 < ^ < 10. In fact, models with ^ < 2 are unstable, while, on 
the other hand, models with ^ > 10 are very similar to each other on the large scale, 
except for the development of a nuclear concentration on scales much smaller than those 
studied in this paper. 

Given the focus of the present paper on the i? 1 /™ photometric models, as a preliminary 
investigation we have compared the photometric profiles of the f x dynamical models 
(treated as "the galaxies" ) to those of the R x l n models (taken as the "fitted models"). 
We have thus obtained the best-fit n(SS?) by minimizing the quantity 



where 77 = R/R c is expressed in units of the effective radius of the models. The 
adopted radial interval is 0.1 < 77 < 10; note that the adopted outer radius is well far out 
with respect to the radius i? ma x typically considered for observed photometric profiles. 

We found that for models with * > 5 a fit over the adopted radial interval can 
be obtained with (A/ii) max < 0.4. In contrast, for \P < 5, in order to obtain a fit with 
(A/i) max < 0.4 we had to reduce the radial interval to typical values 0.23 < 77 < 4. This 
smaller radial range is actually applicable to realistic conditions. In these latter cases 
of low-VP models we found that the best-fit value of n is sensitive to the specific radial 
interval adopted; for example, the quality condition (A/i) max < 0.4 can be satisfied for 
a different radial range (say 0.1 < 77 < 2), but the resulting best-fit value of n turns 
out to be significantly lower. In Fig. A.l we show the optimal fitting relation n = n(ty) 
for the fits based on the full radial range, for \f > 5.5, and on the reduced radial range 
(0.23 < 77 < 4) for ^ < 5.5. Note how close to the de Vaucouleurs profiles, n — 4, are the 
high- 'I' models, as realized from the very beginning (see Bertin & Stiavelli 1993). 




(A.5) 
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Fig. A.l. The best-fit n, obtained by fitting the models, projected along the line 
of sight, with R}l n profiles, as a function of the model dimensionlcss central potential 
\&. The radial interval adopted in the fit is described in the text. Note the "plateau" at 
n = 4 reached by high- 'J models. 

Appendix B: Fitting procedures 

In this Appendix we briefly describe different options to fit a model, chosen within the 
set of models presented in Appendix A, to the photometry of observed objects available 
within a given spatial range. 



B.l. Profiles, circularized profiles, curves of growth 

In the procedure of profile fitting, we generally refer to circularized profiles. The objects 
selected for the present study are fairly round in projection. Therefore, the differences 
between the photometric profiles taken along the short (with coordinate R s ) or the long 
(with coordinate Ri) axis (as preferred by CCD93, DCC94) with respect to the profiles 



obtained by referring to the circularized radius R = ^JR s Ri are rather small. 

The observed profiles in mag/arcsec 2 , n \, s {R) = — 2.5 log 7 bs , are compared to 
the corresponding model profiles J mo d 

MmodOR) = -2.5Iog/ mod (i2). (B.l) 
Following standard methods (e.g., see Bertin, Saglia, & Stiavelli 1988), the best-fit pa- 
rameters are determined by minimizing the % 2 function defined as: 

2 _ [Mobs(-Rj) — Mmod(-Ri)] 2 



X. 



\ "> [Mob 



(B.2) 

o~obs(Ri) represents the 



where N is the number of data-points in the interval (i? m i n , Rm&xj, • 
photometric error on ^ b s (-Ri)> an d the factor (N — N par ) takes into account the number 
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of free parameters in the fit. The formal errors Si associated with the best-fit values of 
the parameters Xi involved can be calculated from the covariance matrix E = 2H~ 1 , 
where H is defined as 



while the off-diagonal terms provide the correlation coefficients. 

The parameters thus determined (e.g., the effective radius R c ) do not represent the 
quantities implied by their ideal definitions (e.g., the projected radius of the disk that 
encloses half of the total luminosity), because the fit is performed only on a limited radial 
range (i? m i n , Rmax)- In this procedure the model parameters (e.g., n, R c , and Iq for the 
R}/ n models) are taken to be independent of each other. 

A second fitting procedure (preferred, e.g., by PS97) refers to the so-called curve 
of growth, i.e. to the function associated with the integrated luminosity L(R) = 
2tt §q R' I(R')dR' . One then compares the data points for m b s (R) = — 2.5 logi b s (-R), 
expressed in magnitudes, to the corresponding curve of growth m mo d(-R) derived from 
the photometric model. It is clear that this procedure involves a major difference with 
respect to the profile fitting procedure, because the data points of the curves of growth 
cannot be treated as independent of each other. If we proceed naively, but improperly, by 
means of a x 2 function defined as in Eq. (B.2), (but in the results shown in Section 2.3 
we took a constant photometric error) with \x replaced by m, we find "best-fit values" of 
the relevant parameters that can be significantly discrepant with respect to those found 
from profile fitting (compare the values of n given in the second column of Table 3 with 
those recorded in Table 1). The discrepancies follow the general trends reported in the 
literature (see Introduction and PS97). 

B.2. Constraints 

One of the tests that we have performed on our set of observed photometries has been 
the study of the dependence of the best-fit parameters on the radial range of the fit. In 
order to compare with earlier studies (e.g., see Burkert 1993), we have thus considered a 
case where all photometric profiles for the four galaxies studied in this paper were limited 
to a common excursion of 4.5 mag (see second group of columns in Table 3). 

Another natural possibility (see Bertin, Saglia, & Stiavelli 1988) is to consider a fit 
under the constraint that the observed profile and the model profile be characterized 
by the same integrated luminosity within the radial range of the fit. To these fits we 



Ha = (N — 7V par ) 




(B.3) 



Therefore, we have 




(B.4) 
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refer as fits with constrained luminosity (see Sect. 2, third group of columns in Table 3). 
Explicitly, we impose: 



respect to which the x 2 function is minimized is reduced by one. 

Similarly, for the fit to the curves of growth it is natural to impose that £ bs(^max) = 



B.3. Comments 

Not unexpectedly, when free parameters are available, such as n or f , the search for 
the best-fit model actually leads to several possible minima of the relevant x 2 function. 
In Table 2 we have recorded the properties (in parentheses) of a secondary minimum 
available for NGC 1379, within the models, which appears to be particularly signif- 
icant. In Tables 1 and 3 we have considered only the R}l n models corresponding to the 
absolute minima found. For each galaxy, we have also checked the properties of other 
(typically 4 or 5) models corresponding to secondary minima, with x 2 within a factor 
of 3-4 with respect to that associated with the absolute minimum. The related model 
parameters turn out to be significantly different only when the associated value of x 2 
changes significantly. The values of n and R e thus identified tend to increase, or decrease, 
together, but the photometric parameter (f>pp remains robustly stable. 

The case illustrated briefly in the second group of columns in Table 3, that of the 
reduced radial range, presents, not surprisingly, the most complex behavior. Here we can 
even witness the occurrence of "crossing" of solutions. An example of this phenomenon 
is provided by the galaxy NGC 4374. For this object the absolute minimum over the 
reduced radial range is obtained for n s=a 14.5 and R c ~ 1092, with x 2 ~ 0.34 (see Table 
3). A very near (when considering the x 2 value) secondary minimum is n m 9.9 and 
R e w 302, with x 2 ~ 0.55. When fitting the profile of the same galaxy over the extended 
radial range the absolute minimum is found (see Table 1) at n w 9.5 with R e w 259 and 
X 2 ~ 0.46. Thus, the solution corresponding to the absolute minimum obtained over the 
the whole available radial range can turn out to be associated with a secondary minimum 
when the fit is extended over a reduced radial range. 

Appendix C: General properties of the velocity dispersion profiles of simple 
dynamical models 

If we wish to draw consequences, from the observed trends in the photometric profiles, 
on the structure of the FP, we have to take into account some dynamical information. 
In practice, for the objects under consideration we have to estimate the relevant virial 



*'-Rmi„ J-Rmin 

so that the luminosity scale is set as a constraint and the number of free parameters with 
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coefficient (see Appendix D). Since the dynamical quantity that appears in the FP is the 
central (projected) velocity dispersion (based on a suitable small aperture), it is clear that 
our consequences will depend mostly on the assumed behavior, for the velocity dispersion 
profiles, in the vicinity of the galactic center. In this Appendix we will show that if we 
start from a given density profile and then proceed to infer, as commonly done (through 
the Jeans equations under simple assumptions on the pressure tensor; e.g., see PS97), 
the velocity dispersion profile, the resulting kinematical profile is very sensitive to the 
specific limiting behavior (at very small radii) of the density profile we consider. In many 
cases the implied kinematical profile may be completely unrealistic (e.g., by developing 
a major "central hole", not observed; see Bailey & MacDonald 1981). This sensitivity is 
unfortunate, because it makes it necessary to address in detail the dynamical basis of the 
photometric models used to fit the observations. In this respect, we should recall here 
that the velocity dispersion profiles of the /oo models (see Appendix A. 3) are (generally 
realistic) monotonic, decreasing with radius. 

In the following we focus on spherically symmetric, isotropic, one-component galaxy 
models with finite mass and, for simplicity, limit our discussion to the unprojected velocity 
dispersion only. 

Suppose that the model density profile at small radii is characterized by the asymp- 
totic behavior 

p(r)~r- a , (C.l) 
with < a < 3. From the explicit solution of the appropriate hydrostatic equilibrium 
condition, 

± WW;ir)] ^ G ^mi , <G2) 
it is readily shown that for a = the velocity dispersion reaches a finite (non vanishing) 
value at the center (dependent on the shape of the entire density profile). For < a < 1 
we find 

a r 2 (r)~r Q , (C.3) 
and so of (0) = and a "central hole" is present (note that all R}l n models fall in this 
case, because close to the center their density profile is characterized by a = 1 — 1/n). 
When a = 1, we have 

a r 2 (r) ~ -rln(r) , (C.4) 
and so the central value of the velocity dispersion vanishes. Finally, for 1 < a < 3 the 
relation is 

cr r 2 (r) ~r 2 - Q . (C.5) 
Thus, for 1 < a < 2 the central velocity dispersion vanishes and a hole occurs; when 
a — 2 the central velocity dispersion is finite (its value depends on the shape of the entire 
density profile, as in the case a = 0), and for 2 < a < 3 it diverges. 
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In order to obtain an estimate of the velocity dispersion gradient close to the center, 
we can rewrite Eq. (C2) as 



Let us now focus on the case a = 0, and thus express the density profile as p(r) ~ 
p(0) — pir 13 . Then, it can be shown that for < (3 < 2 there is indeed a hole in the 
velocity dispersion; in turn, for (i = 2 the direction of the slope of the velocity dispersion 
depends on the density profile of the entire model, while for (3 > 2 the velocity dispersion 
is monotonically decreasing with radius. 

We may add that, independently of the slope of the stellar density profile in the 
central regions, if a massive black hole is present, then the slope of the square of the 
velocity dispersion is -1, and so, at least in the central regions, no hole is present in the 
velocity dispersion profile. 

By experimenting with more general density profiles in the central regions, we have 
also encountered the possibility of "wavy" velocity dispersion profiles. 

Appendix D: The virial coefficient iSTv 

After this short kinematical digression, it should be clear that the virial coefficient Ky = 
GM I (R c <Jq) that controls the relation between the virial theorem and the FP is sensitive 
to the dynamical model that we adopt to describe galaxies, not just to the photometric 
signature as captured, for example, by the R}> n profiles. In addition, within a given set 
of dynamical models, the virial coefficient can be sensitive to the precise definition that 
we adopt for "central velocity dispersion" . 

To set the basis for our simulations of Sect. 3, we have computed the virial coeffi- 
cient for the isotropic i? 1 /" models, with <jq based on an aperture of radius i? c /10. In 
addition, to illustrate the above-mentioned sensitivity problem we have computed the 
same coefficient also for the models. We have checked that, for the isotropic i? 1 /™ 
models, the numbers are in agreement with those reported earlier by PS97. The results 
are summarized in Fig. D.l. 

From this figure, an immediate warning in relation to the sensitivity problem derives 
from inspection of the high- 'I' models: while with respect to the photometric performance 
the models with 7 < W < 10 are all very well fitted by R X I A profiles (see Fig. A.l), their 
virial coefficient is found to vary significantly, from fa 2.3 to w 3.3. 

Another simple way to point out quantitatively the embarassing modeling problem 
that we are facing is given in Table D.l. Here we record, for the four galaxies studied in 
this paper, the value of the virial coefficient obtained by referring to the best-fit i? 1 /™ 
models suggested by Table 1, to those of Table 3 (b: reduced radial range), or to the best- 
fit /oo models identified in Table 2. Note that if we modeled these galaxies as isotropic 




(C.6) 
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Fig. D.l. The virial coefficient for the (triangles) and the isotropic R}l n (squares) 
models based on an aperture of radius Re/ 8. 

i? 1 / 4 models, they would all be assigned the same virial coefficient Ki/ 4 w 4.6. We should 
emphasize that, in the spirit of the studies of the FP, all these inferences are made without 
considering the detailed behavior of the actual velocity dispersion profiles for the galaxies 
under investigation. 

Table D.l. Virial coefficient Ky corresponding to different fitting procedures. 



NGC 


(i? 1/n ) a 


(R 1/n ) b 


foe 


1379 


6.8 


6.8 


3.2 








(5.1) 


4374 


1.8 


1.8 


2.4 


4458 


5.7 


5.5 


3.0 


4552 


1.6 


1.8 


2.5 



a Extended radial range, no L constraint. See Table 1. 
b Reduced radial range, no L constraint. See Table 3. 
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